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Abstract 

Plant chitinases (EC 3.2.1.14) and chitinase-like (CTL) proteins have diverse functions including cell wall biosynthesis and 
disease resistance. We analyzed the expression of 34 chitinase and chitinase-like genes of flax (collectively referred to as 
LusCTLs), belonging to glycoside hydrolase family 19 (GH19). Analysis of the transcript expression patterns of LusCTLs in the 
stem and other tissues identified three transcripts (LusCTU9, LusCTL20, LusCTLll) that were highly enriched in developing 
bast fibers, which form cellulose-rich gelatinous-type cell walls. The same three genes had low relative expression in tissues 
with primary cell walls and in xylem, which forms a xylan type of secondary cell wall. Phylogenetic analysis of the LusCTLs 
identified a flax-specific sub-group that was not represented in any of other genomes queried. To provide further context 
for the gene expression analysis, we also conducted phylogenetic and expression analysis of the cellulose synthase (CESA) 
family genes of flax, and found that expression of secondary wall-type LusCESAs {LusCESA4, LusCESA? and LusCESAS) was 
correlated with the expression of two LusCTLs {LusCTLl, LusCTLl) that were the most highly enriched in xylem. The 
expression of LusCTL19, LusCTL20, and LusCTL21 was not correlated with that of any CESA subgroup. These results defined a 
distinct type of CTLs that may have novel functions specific to the development of the gelatinous (G-type) cellulosic walls. 



•0-PLOS I o-^E 



Citation: Mokshina N, Gorshkova T, Deyholos MK (2014) Chitinase-LIke {CTL) and Cellulose Synthase {CESA) Gene Expression in Gelatinous-Type Cellulosic Walls of 
Flax {Linum usitatissimum L.) Bast Fibers. PLoS ONE 9(6): e97949. doi:10.1371/journal.pone.0097949 

Editor: Olga A. Zabotina, Iowa State University, United States of America 

Received January 23, 2014; Accepted April 26, 2014; PublisKied June 11, 2014 

Copyriglit: © 2014 Mokshina et al. This Is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction In any medium, provided the original author and source are credited. 

Funding: This work was partially supported by the Russian Foundation for Basic Research (project no. 12-04-31418) and a grant of the President of the Russian 
Federation "Support of Leading Scientific Schools (NSh_825.2012.4)" and by Genome Canada and the Natural Sciences and Engineering Research Council 
(Canada). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 

* E-mail: deyholos@ualberta.ca 



Introduction 

Flax [Linum usitatissimum L.) phloem fibers are a valuable 
industrial feedstock and are also a convenient model system for 
studying secondary cell wall formation. The mechanical properties 
of bast fibers are largely dependent on the composition of their 
secondary walls. Bast fibers have gelatinous-type walls, which are 
rich in cellulose (up to 90%) and lack detectable xylan and ligniu. 
Gelatinous fibers are widespread in various land plant taxa, but 
have been studied primarily in angiosperms. Depending on the 
species, either phloem or xylem (of either primary or secondary 
origin) can produce gelatinous fibers in various organs including 
stems, roots, tendrils, vines, and peduncles [1,2]. The mechanisms 
of gelatinous cell wall development in these fibers remain largely 
unclear. However, some genes imphcated in gelatinous cell wall 
development have been identified. The list includes fasciclin-like 
arabinogalactan proteins (FLAs) [3-6], fi-galactosidases [7,8], and 
lipid transfer proteins [6]. A role for P-galactosidases in G-type 
wall development has been demonstrated functionally [8]. 

Transcripts of genes encoding chitinase-like proteins are 
reportedly enriched in fibers, particularly during the cell wall 
thickening stage of flax phloem cellulosic fiber development [6] . 
Expression of CTLs during primary or secondary cell wall 
deposition has also been reported in species other than flax 
[9, 1 0] . Plant chitinase-like proteins have been identified in a wide 



range of organelles and tissues, including the apoplast and vacuole 
[11]-^^ 

Chitinase-like proteins belong to a large gene family that 
includes genuine chitinases (i.e. proteins with proven chitinase 
activity) and other homologous proteins, which may not have 
chitinase activity [12-15]. Here, we refer to both genuine 
chitinases and their homologs collectively as chitinase-like proteins 
(CTLs). 

Chitinases catalyze cleavage of P-l,4-glycoside bonds of chitin 
and are organized in five classes (Classes I-V), which can be 
distinguished on the basis of sequence similarity [11,16,17]. 
Classes I, II, and IV belong to glycoside hydrolase family 19 
(GH19), found primarily in plants, whereas Classes III and V 
belong to glycoside hydrolase family 1 8 (GH 1 8) present in various 
types of organisms [18-20]. The Class I chitinases are found in 
both monocots and dicots, while classes II and IV are found 
mainly in dicots [21]. Class I and IV chitinases contain a highly- 
conserved cysteine-rich domain - the chitin binding domain 
(CBD) - at the N- terminal region [21,22], but there are two 
characteristic deletions in the main catalytic domain in Class IV 
chitinases [21]. Because chitin is the major component of fungal 
cell walls, chitinases are classic pathogeiiesis-related proteins 
involved in non-host-specific defense [23,24]. 

Plants also contain chitinase-like proteins that are not induced 
by pathogens or stresses. In many cases, these chitinase-like 
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Table 1. Description of plant samples. 





Sample 


Basic tissues 


The main difference from other tissues 


"Apex" 


epidermis, meristem cells, parenchyma cells, procambium, sieve elements, companion cells, 
fast growing phloem fibers with primary cell wall, xylem vessels with primary cell wall and 
secondary thickness appears 


phloem fibers with primary cell wall 


"TOP" 


epidermis, parenchyma cells, sieve elements, companion cells, fast growing phloem fibers with 
primary cell wall, xylem fibers and vessels secondary thickness (xylan cell wall) 


phloem fibers with pnmary cell wall 


"MID" 


epidermis, parenchyma cells, sieve elements, companion cells, cambium, phloem fibers with 
gelatinous cell wall, xylem fibers and vessels with xylan cell wall 


phloem fibers at the beginning of gelatinous 
cell wall deposition 


"Peel" 


epidermis, parenchyma cells, sieve elements, companion cells, fibers with thick gelatinous cell wall 


phloem fibers with thick gelatinous cell wall 


"Xylem" 


parenchyma cells, cambium, xylem fibers and vessels with xylan cell wall 


xylem fibers and vessels with xylan type of 
cell wall 


"Fiber" 


isolated fibers with thick gelatinous cell wall 


isolated phloem fibers with thick gelatinous 
cell wall 


"Root" 


epidermis, parenchyma cells, sieve elements, companion cells, xylem fibers and vessels with xylan 
cell wall 


xylem fibers and vessels with xylan type of 
cell wall 



doi:10.1371/journal.pone.0097949.t001 



proteins have been shown to lack detectable chitinase activity. 
Chitinase-like proteins may play an important role during normal 
plant growth and development [13,1 5,25] . For example, AtCTLl is 
constitutively expressed in many organs oiArabidopsis. Mutations of 
AtCTLl lead to ectopic deposition of lignin in the secondary cell 
wall, reduction of root and hypocotyl lengths, and increased 
numbers of root hairs [15]. It was suggested that this gene could be 
involved in root expansion, cellulose biosynthesis, and responses to 
several environmental stimuli [13,26,27]. In particular, co- 
expression of some CTLs with secondary cell wall cellulose 
synthases (CESAs) was reported [28]. It has been suggested that 
these chitinase-like proteins could take part in ceUulose biosyn- 
thesis and play a key role in establishing interactions between 
ceUulose microfibrils and hemiceUuloses [14]. 

The xylan-type secondary wall is the most common secondary 
wall in land plants and is characteristically rich in cellulose, xylan, 
and ligTiin [2]. Compared to typical xylan-type secondary walls, 
gelatinous layers are enriched in ceUulose, have a higher degree of 
ceUulose crystallinity, larger crystaUites, and a distinctive set of 
matrbc polysaccharides (see [2] and references therein). Presum- 
ably, ceUulose synthase genes have a significant role in gelatinous 
ceU waU formation, but the expression patterns of the complete 
flax CESA famUy has not been described to date. It is known that at 
least three isoforms of CESAs comprise the cellulose synthase 
rosette: CESAl, CESAS, and CESA6 are required for cellulose 
biosynthesis in primary ceU walls [29], whereas CESA4, CESA7, 
and CESAS are required for cellulose biosynthesis during 
secondary waU deposition [30] . 

Flax is a useful model for comparative studies of cell waU 
development: different parts of the flax stem form a primary ceU 
wall, xylan type secondary ceU waU, or gelatinous ceU waU; these 
stem parts may be separated and analyzed by diverse approaches, 
including functional genomics. Furthermore, the publication of a 
flax whole genome assembly [31] facilitates a thorough study of 
key gene famihes. 

In the present study, we measured expression of all predicted 
LusCTL genes of the GH19 famUy in various tissues including 
those that produce gelatinous-type and xylan-type ceU waUs. We 
also described the LusCESA gene famUy and measured expression 
of its transcripts in comparison to LusCTLs. Phylogenetic analysis 
of Lu sCTL and LusCESA genes identified distinct groups oi LusCTL 
genes that were expressed preferentially at specific stages of bast 
fiber gelatinous ceU waU development. 



Materials and Methods 

Plant Growth 

Flax {Linum usitatissimum L.) var. Mogilevsky plants were grown in 
pots in a growth chamber at 22°C, with a light intensity of 
approximately 200 |lE on a 16 h light/8 h dark cycle. Plants were 
harvested at the period of rapid growth (4 weeks after sowing). Plant 



Apical part 
of stem 




Figure 1. Scheme of sample collection. The segments of the flax 
stem (apical part "Apex", top region TOP, middle region IVIID, Peel 
containing phloem fibers, Xylem and isolated Fibers) were taken for 
RNA extraction. A detailed description of each sample is presented in 
Table 1. Bar: 50 |im for "Apex", 100 |im for TOP and MID, 1 mm for 
Fibers. 

doi:10.1371/journal.pone.0097949.g001 
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material was sampled with respect to the location of the snap point, 
which is a mechanically defined stem position in which fibers 
undergo transition fi-om elongation to secondary cell wall formation 
[32] . The following seven samples were collected (Table 1 , Figure 1 ): 
1. "Apex" - the apical part of stem (1 cm of length). 2. "TOP" - the 
following "apex" segment of stem above the snap point with phloem 
fibers in the process of elongation. 3. "MID" - the stem segment 
(5 cm of length) below the snap point which contained fibers at early 
stages of secondary cell wall thickening. 10 cm of the stem 
downwards from "MID" was divided into Peel (4), which contained 
epidermis, parenchyma cells, phloem fiber bundles and sieve 
elements and Xylem (5), which contained parenchyma cells, xylem 
vessels and xylem fibers. 6. "Fibers" - i.e. isolated phloem fibers 
were obtained by washing Peels in 80% ethanol in a mortar several 
times and gently pressing the fiber-bearing tissues with a pestie to 
release the fibers. 7. Roots. The number of biological replicates was 
three, with five plants in each replicate. 

Sequence Alignment and Phylogenetic Analysis 

Predicted amino acid and nucleotide sequences of CTLs (Pfam 
domam: PF00182) and CESAs (PF03552) were obtained from the 
Phytozome database v.9.0 (Limm usitatissimum, Populus trichocarpa, 
Arabidopsis thaliana). CESAs of poplar (PtiCESAs) were renamed 
according to Kumar et al. [33]. A list of various well-characterized 
CTLs from different plant species was obtained from previously 
published works [10,21]. Sequences were ahgned using MUSCLE 
with default parameters, and a phylogenetic tree was constructed 
using MEGA5 based on the Maximum Likelihood and Neighbor- 
Joining methods [34], bootstrapping 1000 replicates [35], model 
WAG+G or JTT+G. Signal peptides for protein sequences were 
predicted using SignalP (http://www.cbs.dtu.dk/services/SignalP/ 
), molecular weights, isoelectric points of the proteins were analyzed 
by ProtParam (http:/ /web. expasy.org/protparam/). 

Reverse Transcription Quantitative Real Time PCR 

Total RNA from all plant samples was isolated using a Trizol- 
extraction method combined with an RNeasy Plant Mini Kit 
(Qiagen) according to the manufacturer's instructions. RNA 
quality was evaluated by electrophoresis using a BioAnalyzer 
(Agilent), and no degradation of RNA was evident. Residual DNA 
was eliminated by treatment with DNAse I using the DNA-free kit 
(Ambion). Gene specific primers for CTL and CesA genes were 
designed using Universal ProbeLibrary Assay Design Center 
(http://www.roche-applied-science.com/shop/CategoryDisplay? 
catalogid = lOOOl&tab — &identifier = Universal+Probe+Library) 
(File SI). One microgram of total RNA was reverse-transcribed 
using RevertAid H Minus First Strand cDNA Synthesis Kit 
(Thermo Scientific). The cDNAs were diluted 1:32 with nuclease 
free water. Real-time PCR was performed in a 7900 HT Fast real- 
time PCR system (Applied Biosystems, USA). Each 10 |J.L real- 
time PCR cocktail contained 2.5 |J,L of 0.4 |J.M concentrations of 
both forward and reverse gene-specific primers, and 2.5 |J,L of 
cDNA, 5 |J.L of 2 xDynamite qPCR mastermix (Molecular Biology 
Service Unit - University of Alberta) which included SYBR green 
(Molecular Probes) and Platinum Taq (Invitrogen). The thermal 
cycling conditions were 95 °C for 5 minutes, 40 cycles of 95°C for 
15 seconds, and 60°C for 1 minute. A 60-95°C melting curve was 
performed to confirm the specificity of the products. Threshold 
cycles (CT) were determined using 7900 Fast Software. Cx values 
were normalized using eukaryotic translation initiation factors lA, 
5A [LusETIFl , LusETIFSA) and glyceraldehyde 3-phosphate 
dehydrogenase (LusGAPDH) gene from flax (File SI) [36]. From 
each of three biologically independent cDNA samples, two 
independent technical replications were performed and averaged 
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Figure 2. Phylogenetic relationships among LusCTLs. All predicted LusCTLs amino acid sequences were aligned using their deduced full-length 
peptide sequences In the MEGA platform (MUSCLE), Maximum Likelihood Method, JTT+G model, bootstrap 1000. 
doi:1 0.1 371/journal.pone.0097949.g002 



for further calculations. AACT values were generated using the 
apex sample as a reference. Relative transcript abundance 
calculations were performed using comparative Ct (ACt) method 
as previously described [37] for flax tissues (TOP/Apex, MID/ 
Apex, Peel/ Apex, Xylem/Apex, Fiber/ Apex, Root/ Apex). Heat 
maps of expression levels of some genes were then created with 
MeV v4.8 (Multi Experiment Viewer, http://www.tm4.org/mev.) 
using the means of ACt. 

Results 

LusCTL Phylogenetic Characterization 

We searched within the flax genome assembly (version 1.0) for 
predicted genes with homology to Pfam domain PF00182, which is 



characteristic of chitinases of the glycosyl hydrolase family 1 9 (GH 
19) family [31,38]. This search identified 37 predicted chitinase or 
chitinase-like genes (referred to here collectively as LusCTLs) 
(Table 2). However, only three of the predicted proteins {LusCTLS, 
LusCTLI 0, and LusCTLI 5) contained a conserved chitin-binding 
domain (CBD), suggesting that not all of the LusCTLs use 
chitinase as a substrate. The mean predicted protein size of the 37 
LusCTLs was 246.5 aa (or 27 kDa), and the majority (30/37) were 
predicted to be secreted (Table 2). 

The labels assigned to the 37 predicted LusCTLs are shown in 
Table 2. LusCTLI and LusCTL2 were so named because they 
encoded proteins that were most similar to CTLl and CTL2, 
respectively, which have been characterized in other species (e.g. 
A. thaliana [14] and G. hirsutum [10]) (Table 2). The gene LusCTL37 
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Figure 3. Phylogenetic relationship among LusCTLs and selected CTLs of 43 different plant species. All the predicted CTLs sequences 
were aligned using their deduced full-length peptide sequences in MEGA (MUSCLE). Maximum Likelihood Method, WAG+G model, bootstrap 1000. 
Classes l-IV were labeled according to Hamel et al. [21] and Zhang et al. [10] and Groups A-C are as shown in Figure 2. 
doi:1 0.1 371/journal.pone.0097949.g003 
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Figure 4. Relative differential expression of LusCTL genes in different tissues of the flax stem. AACj values for each of CTLs were 
generated using an average ACj for all CTLs as a reference. Error bars show the standard error of the mean. 
doi:1 0.1 371/journal.pone.0097949.g004 
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Figure 5. Phylogenetic relationship among LusCESAs, AtCESAs, 

PtiCESA. All CESAs sequences were aligned using their deduced full- 
length peptide sequences on MEGA platform (MUSCLE). Maximum 
Likelihood Method, JTT+G model, bootstrap 1000. 
doi:1 0.1 371 /journal.pone.0097949.g005 

was predicted to encode only a protein of 69 aa, whicli is much 
shorter than the rest of the LusCTLs (Table 2), and so it was not 
used in further analyses. 

The LusCTLs and their inferred phylogenetic relationships are 
shown in Figure 2. Based on this dendrogram, the predicted 
LusCTLs were divided into three groups: Group A included 
LusCTLl-6, Group B included LusCTL7-14, and Group C 
included LusCTL15-36 (Figure 2, Table 2). The flax-specific tree 
shown in Figure 2 was expanded by the addition of representative 
GH19 CTLs from other species (Figure 3). In this multispecies 
tree, LusCTLs of Group A, which includes LusCTLl and 
LmCTL2, were part of the same clade as the well-characterized 
AtCTL2 of A. thaliana and GhCTLl, GhCTL2, GhCTLVII of G. 
hirsutum. The Group B LusCTLs {LusCTL7-14) were in the same 



clade as the previously defined Classes I, II, III, GH19 chitinases 
[10,21]. Most of group B was in the same sub-clade as Class II, 
although none of the previously defined Classes I-III were 
monophyletic in our analysis. Finally, our Group C LusCTLs 
[LusCTLl 5-36] formed a monophyletic clade with representatives 
of the previously defined Class IV GH19 chitinases. 

LusCTL Transcript Expression 

Quantitative real-time reverse-transcription PCR (qRT-PCR) 
was performed to study LusCTL expression patterns of L. 
usitatissimum genes of chitinase-like proteins in various tissues and 
stages of development (Table 1). These tissues and their names as 
used here are equivalent to the names used in previous studies 
[6,32,39,40]. Only 34 sets of primers were used in this assay, 
because members of each of two pairs of LuCTLs could not be 
distinguished by unique primers: LusCTL28 and LusCTL29 (95.6% 
aa and 96.3% nt identity), and LusCTL33 and LusCTL34 (99.6% 
aa and 98.8% nt identity). Thus a common set of primers was used 
for each of these pairs. We observed that transcripts of LusCTLl 
showed enriched levels of expression (compared to the apical part 
of stem) in tissues in which cell walls were undergoing thickening 
(Figure 4) in xylem and in phloem fibers. Transcripts for this gene 
were enriched .57-fold in xylem, 28-fold in the MID region, and 
20-fold in fiber. Another predicted CTL, LusCTL2, showed a 
similar pattern of enrichment in secondary-wall bearing tissues 
(8.3, 4.5 and 1.4-fold higher in xylem, MID and fiber, respectively, 
compared to the apex), although the magnitude of its enrichment 
was not as strong as LusCTLl. These two LusCTLs had high 
sequence similarity to each other (91.9% amino acid identity) and 
had similar patterns of expression as compared to each other in the 
various flax tissues. 

A subset of LusCTL genes [LusCTLl 0, LusCTLl 1, LusCTLl 9, 
LusCTL20, LusCTL21, LusCTL23, LusCTL24, LusCTL2G) had high 
relative expression in tissues that contained phloem fibers (MID, 
peel, fiber) but low relative expression in xylem (Figure 4). Three 
of these genes [LusCTLl 9, LusCTL20, LusCTL21) were enriched > 
40-fold in fibers compared to the apical part of stem. These three 
genes had high similarity to one another (76% identity between 
LusCTL19 and LusCTL20 as well as between LusCTL19 and 
LusCTL2T, 91% identity between LusCTL20 and LusCTL21) 
(Figure 2). 

LusCESA Phylogenetic and Expression Characterization 

To provide context for the expression patterns of LusCTLs, and 
to test whether the expression pattern of cellulose synthase 
(LusCESA) genes differed between gelatinous fibers and cells with 
a xylan type of secondary cell wall, expression of LusCESAs in 
difierent flax tissues was analyzed. We identified 14 predicted 
LusCESAs in the flax whole genome assembly by searching 
predicted proteins for the conserved cellulose synthase domain 
(Pfam PF03552). No putative LusCESA? genes were found in the 
original published genome published (vl.O, [31]). However, 
though BLAST alignment of the CDS of Arabidopsis and poplar 
CESA7 sequences, two scaffolds (scafibld_57 and scafibld_464) of 
the flax genome assembly were identified as encoding CESA7 
homologs, and these were annotated using the Augustus server 
(http://bioinfuni-greifswald.de/augustus/). Thus, all 16 predicted 
LusCESAs were aligned with well-characterized AtCESAs from A. 
thaliana and PtiCESAs from P. trichocarpa (Figure 5). This alignment 
was used to construct a phylogenetic tree and annotate the 
LusCESAs, which were named according to the established A. 
thaliana [41] and P. trichocarpa nomenclature systems (Table 3, 
[33]). The number oi LusCESAs and PtiCESAs isoforms identified 
for each of the eight major types of CESAs was similar except in 
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Table 3. Predicted cellulose synthase sequences In Linum usitatissimum. 





Locus I.D. 


Code 


Pti homolog 


At homolog 


Length, aa 


MW, kDa 


Pl 


Lusl0018902 


LusCESAl-A 


PtiCESAl-A, 


AtCESAl 


1079 


121.7 


6.3 


Lusi 0028597 


LusCESAI-B 


PtiCESAI-B 


(AtCESAlO) 


1079 


121.7 


6.3 


Lusi 0039607 


LusCESA3-A 


PtiCESA3-A, 


AtCESA3 


1069 


120.2 


8.0 


Lusi 0007538 


LusCESA3-B 


PtiCES3-B, 




1092 


122.3 


7.5 


Lusl0012198 


LusCESA3-C 


PtiCES3A-C, PtiCES3A-D 




1094 


122.6 


7.0 


Lusi 0008225 


LusCESA4 


PtiCESA4 


AtCESA4 


987 


111.3 


6.9 


Lusi 00061 61 


LusCESA6-A 


PtiCESA6-A, 


AtCESA6 


1074 


121.4 


6.8 


Lusi 0041 063 


LusCESA6-B 


PtiCESA6-B, 


(AtCESA2, 


1096 


123.7 


6.8 


Lusi 0003526 


LusCESA6-B 


PtiCESA6-C, 


AtCESAS, 


1097 


123.7 


6.8 


Lusi 0002939 


LusCESA6-D 


PtiCESA6-D, 


At CESA9) 


1097 


123.5 


7.1 


Lusi 0002940 


LusCESA6-E 


PtiCESA6-E, 




906 


102.1 


7.3 


Lusi 0022449 


LusCESA6-F 


PtiCESA6-F 




1035 


116.9 


6.1 


scaffold_57 


LusCESA7-A 


PtiCESA7-A, 


AtCESA7 


978 


111.2 


6.8 


scaffold_464 


LusCESA7-B 


PtiCESA7-B 




959 


108.3 


8.6 


Lusi 0007296 


LusCESA8-A 


PtiCESA8-A, 


AtCESAS 


988 


111.4 


6.4 


Lusi 0029245 


LusCESA8-B 


PtiCESA8-B 




988 


111.3 


6.7 



doi:10.1371/journal.pone.0097949.t003 



the case of CESA3, where one more gene was identified in P. 
trichocarpa than in L. usitatissimum (Figure 5, Table 3). The LtisCESA 
appeared to be typical of other genes in this family in that they 
were large integral membrane proteins with eight predicted 
transmembrane domains, a hydrophilic domain that faces the 
cytosol, and a zinc fmger domain at the N-terminus of proteins 
with the characteristic amino acid motif "CxxC" (specific for 
CESAs only [42]). 

Relative differential expression of LusCESA genes in different 
tissues of the flax stem was estimated (Figure 6). LusCESA4, 
LusCESA8~A, LusCESAS-B, LusCESA7-A, LusCESA7-B had high 
expression in tissue that produce secondary walls (TOP, MID, 
Xylem, Fiber, Root). Transcripts of these LusCESA isoforms were 
the most enriched in Xylem, which contained cells with xylan-type 
cell walls, and in roots, where secondary vascular tissue (xylem) 
was also well-developed. These secondary cell wall type LusCESAs 
had also high relative expression in ceUulosic fibers, although it 
was not as strong as for xylem. 



Changes in expression of the LusCESA4, 7, 8 isoforms and 
"xylem-specific" LusCTLl and LusCTL2 were well-correlated in 
different flax tissues (Figure 7 A). This group of genes was highly 
expressed in tissues with secondary cell walls (MID, Xylem and 
Roots). In contrast, the "fiber-specific" LusCTLs had very different 
patterns of expression in the same tissues (Figure 7B): these had 
low level of expression in xylem, but high level of relative 
expression in tissues with gelatinous fibers (peel and fiber). 

Discussion 

Certain fibers of many plant species form G-type cell walls, 
which are rich in crystalline cellulose [1]. Expression of CTLs has 
been previously reported to be enriched during development of G- 
type cell walls, along with specific FLAs [4,6,7], LTPs [6] and 
BGALs [6—8] . In this work, we analyzed expression of all LusCTL 
genes of GH 19 in different flax tissues and compared this 
expression with LusCESAs and to their inferred phylogenies. 




>^ y ^''^ 



& V ^ ^ ry <f' ^ ^ ^ 



V 



V 



V V 



■ TOP/Apei 

■ MID/Apex 

■ Peel/Apex 

■ Xylem/Apcx 

■ Fiber/Apex 
Root/Apex 



Figure 6. Relative expression of LusCESAs genes in different tissues of the flax stem. AACj values were generated using the "Apex" sample 
as a reference. Error bars show the standard error of the mean. 
doi:1 0.1 371/journal.pone.0097949.g006 
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LusCESA4 I 
LusCESA7-A ^ 
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LusCESA7-B l 
LusCESS-A 
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Figure 7. Expression of secondary wall LusCESAs and LusCTLs. A - Relative differential expression of secondary cell wall LusCESAs 4, 7A, 7B, 
8A, 8B and LusCTLl and LusCTL2 in different tissues of the flax stem and centroid graph (from the left) for the mean of ACj of these 7 genes. B - 
Relative differential expression of fiber-enriched LusCTLs {LusCTLlO, 11, 19, 20, 21) in different tissues of the flax stem and centroid graph (from the 
left) for the mean of AQ of these 5 genes. The mean of ACj of LusCTLs and LusCESAs were used to produce a heat map. 
doi:10.1371/journal.pone.0097949.g007 



In the flax genome, 16 predicted LusCESAs were identified 
(Table 3). Previously only partial sequences of some flax CESAs 
were published [43]. All 16 flax CESAs could be placed in discrete 
clades with Arabidopsis and Populus CESA homologs (Figure 5). We 
generally numbered LusCESAs in a way that reflects the association 
of each flax gene with its nearest relative in the Arabidopsis genome, 
as was done for CESAs of Populus [33]. Following this pattern, the 
LusCESA6A—F genes we named as a group, similar to PtiCESA6A—F 
and were not distinguished as CESA2/ 9/5/6 as in Arabidopsis clade 
(Table 3) [44] . Most of the flax and Populus CESA genes are present 
as pairs of paralogs in their respective genomes, although there 
were three LusCESAS genes {LusCESA3A-C) for only two Populus 
genes and one Arabidopsis gene. AtCESAl and AtCESAlO were 
represented by only one pair of genes {LusCESAlA, B) in flax. 

It is well established that proteins encoded by different sets of 
three CESA genes (CESAl, 3, 6 and CESA4, 7, 8) are required for 
cellulose synthesis during primary and secondary wall formation, 
respectively [44—46] . The functional relationships of the various 
paralogs of LusCESAs (except LusCESA4, Table 3) are presentiy 
unclear. According to the data obtained here, secondary cell wall 
LusCESA4, LusCESA7-A, B and LusCESAS-A, B were highly 
expressed both in the xylem cells with lignified cell walls (i.e. 
xylan type) and in the phloem fibers with thick gelatinous cell wall. 
This suggests that phloem fibers and xylem may use similar, rather 
than specialized rosettes. This is consistent with observations from 
poplar showing only minor differences in expression of cellulose 
biosynthetic genes in tension wood as compared to normal wood 
[3] . The different properties of gelatinous and xylan type cell walls 
are therefore likely determined not by CESAs, but by other 
proteins associated with cellulose synthesis, which could include 
specific CTLs. 

We observed two LusCTLs that were expressed more strongly in 
xylem tissue than in any other tissue surveyed (Figure 3, LusCTLl, 
LusCTLl). The co-expression of certain isoforms of LusCTLl, 
LusCTL2 and the secondary wall LusCESAs {CESA4, 7, 8) suggested 
a role for these LusCTLs in secondary cell wall development 
(Figure 7). As noted above, LusCTLl and 2 are highly similar to 
AtCTL2 o{ A. thaliana and GhCTLl, GhCTL2, of G. hirsutum. The 
role of CTL2, and its close homolog CTLl, in cell wall biosynthesis 



is especially intriguing since associations between CTLs and 
primary or secondary cell waU synthesis have been reported in 
different plant species [10]. CTL2 is strongly upregulated during 
secondary wall formation in interfascicular fibers in A. thaliana. 
Reduction in crystalline cellulose content in ctll ctl2 mutants was 
demonstrated, leading to the to the suggestion that AtCTLs are 
involved in cellulose assembly. Furthermore, in P. trichocarpa, 
expression of chitinase genes related to AtCTLl, AtCTL2, and 
GhCTLVII are highly correlated with secondary wall formation of 
xylem [47]. It has therefore been proposed that CTLl and CTL2 
work in conjunction with primary- and secondary-cell wall 
CESAs, respectively [14]. One of the hypotheses for CTLl/2 
function is regulation of cellulose assembly and of interaction with 
hemicelluloses via binding to emerging cellulose microfibrils [14]. 
However, the mechanism of CTL action in cell wall biosynthesis 
as well as substrates of catalytic activity (if any) remains unknown. 
It was suggested that the likely substrates of plant chitinases may 
be arabinogalactan proteins, chitooligosaccharides and other 
GlcNAc-containing glycoproteins or glycolipids [13,15,48,49] 
and the mechanism by which CTLs act is more likely to involve 
binding of chitin oligosaccharides than catalysis. Also, it is assumed 
that chitinases may participate in the generation of such signal 
molecules that regulate the organogenesis process [50]. 

Although relative expression oi LusCESA [4, 7, 8) and LusCTLl, 
LusCTL2 in xylem tissue was higher compared with phloem fibers, 
we cannot exclude involvement of these LusCTLs in phloem fiber 
cell wall development. At the same time, a distinct group of 
LusCTLs [LusCTLl 9, LusCTL20, LusCTL21) had very high 
enrichment in samples with phloem fibers (MID, peel, fiber) with 
a low level of expression in xylem. According to the phylogenetic 
tree, these LusCTLs (group C) were most similar to the previously 
defined Class IV chitinases (Figure 4). High constitutive expression 
of Class IV (along with Class I) in most organs of ^. thaliana under 
normal growth conditions has been previously noted [51]. 
Detailed bioinformatic characterization of genes of LusCTL 
distinct group should be conducted in future. Probably LusCTLs 
that are highly expressed in fibers may be specific to the gelatinous 
cell wall, while LusCTLl and LusCTL2 are essential for wall 
thickening in general. 
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Conclusion 

High expression of specific LusCTLs was obsen^ed in diffcix-nt 
types of thick cell wall producing tissues. LusCTLl and LusCTL2 
were preferentially expressed during secondary wall deposition of 
xylem and were coexpressed with secondary cell wall CESAs [4, 7, 
8). Another group of LusCTLs, (especially LusCTLl 9, LusCTL20, 
LusCTL21) were highly expressed in bast fibers, which have 
ceUulose-rich, gelatinous walls. The group of fiber-enriched 
LusCTLs was expanded in flax compared to species that do not 
produce bast fibers, suggesting that these genes might play a 
unique role during gelatinous cell wall development in general and 
cellulose synthesis in particular. It is possible that the presence of 
fiber-specific LusCTLs, along with other key participants, deter- 
mines differences in mechanisms of xylan and gelatinous cell wall 
formation. To estabhsh the functions of these LusCTLs further 
characterization, including analysis of enzyme activity and 
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